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Abstract 

We describe the behavior of two coupled Bose-Einstein condensates in time-dependent (TD) 
trap potentials and TD Rabi (or tunneling) frequency, using the two-mode approach. Starting 
from Bloch states, we succeed to get analytical solutions for the TD Schrodinger equation and 
present a detailed analysis of the relative and geometric phases acquired by the wave function of 
the condensates, as well as their population imbalance. We also establish a connection between the 
geometric phases and constants of motion which characterize the dynamic of the system. Besides 
analyzing the affects of temporality on condensates that differs by hyperfine degrees of freedom 
(internal Josephson effect), we also do present a brief discussion of a one specie condensate in a 
double- well potential (external Josephson effect). 
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I. INTRODUCTION 



In recent years, concepts which has been restricted to foundation of quantum mechanics 
have been considerably enlarged by spreading out to different domains of physics. With 
the introduction of measurements such as separability 0,0, Q and concurrence , and 
the wider understanding that entanglement, and so nonlocality, is at the core of many-body 
phenomena as quantum-phase transition [fjj], superconductivity [7| and Bose-Einstein con- 
densation 8j, we are considerably far from the time when entanglement and nonlocality were 
confined to fundamental aspects of quantum mechanics. On the other hand, the experimen- 
tal techniques developed over the last decades for manipulating atom-field interaction have 
enabled the building of macroscopic atomic ensembles and the experimental verification of 



fundamental concepts in macroscopic scales 



9|. It is worth mentioning the rapid growth of 



quantum information theory which has conferred to its basic ingredients — the phenomena 
of superposition of states and decoherence, entanglements and nonlocality — a great deal 
of advances towards the accomplishments of quantum logical devices. 

Among the standard tools to generate and detect multipartite entanglements, experi- 
ments in Bose-Einstein condensates (BECs) in dilute gases have deepened our incursion 
towards the quantum nature of macroscopic systems. In particular, experiments with a 
trapped gas of 87 Rb atoms with two different hyperfine sublevels prompt the engineering 
of a Josephson-like coupling between two condensates by a laser-induced Raman transi- 



tion 
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Such "internal" Josephson effect |12j mediates intraspecies collisions apart 



from interspecies one. These atom-atom interactions empowered the investigation of the 



dynamics of the relative phase of coupled condensates 13] and Rabi oscillations [14| on 
macroscopic s yste ms. Moreover, precise measurements of scattering lengths has also been 
accomplished 15|, aiming to quantify properly the non- linear dynamics associated with col- 
lisions. In the two-mode approximation, the coupled condensates has been employed to 
investigate entanglement dynamics (si, 16. \y\ 
detect macroscopic superposition states 



and the possibility to prepare, control and 
" jl8 . 3, 0|. Beyond this achievements, the analysis 



of the Josephson effect in this two-mode exactly soluble model may provide a clue for the 
examination of macroscopic coupling arising in less tractable form of the general theory of 



BECs 



12|. 



Whereas in real experiments the trap potential may be considered to be a time- 
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independent function, excepting for small fluctuations, time-varying scatterin g le ngths are 



usually produced through Feshbach resonances while, as pointed out in Ref. [12J, the am- 
plitude and phase of the laser field may vary in time. In this connection, the present paper 
is devoted to the TD version of the two- mode Hamiltonian (TMH), where the effective fre- 
quency of the trap potential for both atomic species are TD functions, as well as the Rabi 
frequency and the scattering lengths. A similar approach was employed in Refs. 2ll . 
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where, however, only the phase of the external field inducing the Raman transition was 
assumed to be a TD slowly varying function. Instead, our treatment considers time depen- 
dence of all Hamiltonian parameters, focusing on two particular subjects: the analyzes of 
geometric phase acquired by wave function of the whole system and the control of the dy- 
namics of pseudo-spin states governed by the TD TMH. Starting from Bloch states, whose 
preparation is achieved by applying a laser pulse to atoms condensed in a single hyperfme 
level jis| . we demonstrate that its evolution, visualized as a vector on the Bloch sphere, can 
be used to control the geometric phase and the population imbalance following from the 
whole wave function of the condensates. Our treatment also permits a detailed analyzes of 
the relative phase between the condensed states. 



In Ref. 



231 ]. the authors studied the dynamics of a strongly driven two-coupled BECs 



in two spatially localized modes of a double-well potential, where the tunneling coupling 
between the two modes is periodically modulated. In our work we also study the TD TMH 
associated to the "external" Josephson effect, analyzing its differences with relation of the 
"internal" Josephson effect. 

Similarly to the above mentioned fundamental phenomena, the geometric phase has over- 
taken its striking rule on fundamental physics to widening our understanding of phenomena 



as quantum Hall effect 24], |25|, |26( and for the implementation of fault-tolerant quantum 
gates [27(]. After its discovery by Berry on adiabatic processes [281 ] . it has been generalized 



to nonadiabatic [29J, noncyclic [30( and nonunitary [3jJ, |32j quantum evolutions. Recently, it 
has been investigated m^ different areas of physics, ranging from BECs 21] and cavity quan- 



tum electrodynamics 
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341 ] to condensed matter [35|] and quantum information theory [27j . 



In particular, the Berry phase of mesoscopic spin in Bose-Einstein condensates, induced by 
a TD slowly varying driven field, has been investigated under the TMH 2l|, |22J. In our 
treatment, the evolution of the geometric phase of this mesoscopic system is evaluated in a 
more general scenario, where all the Hamiltonian parameters are assumed to be TD. 
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The paper is organized as follows. In Sec. II we introduce and solve the Schrodinger 
equation associated to the TD TMH, presenting the evolution operator. The dynamics 
of BECs for initial Bloch states is analyzed in Sec. Ill, where we show that they remain 
as Bloch states apart from a global phase factor accounting for the elastic collision terms. 
The geometric phase acquired by the state vector of the system is presented in Sec. IV 
and a detailed analyzes of its time evolution is found in Sec. V for different regimes of 
the parameters. In Sec. VI we take the problem of a TMH from a different perspective, 
considering the external Josephson effect instead of the internal one. Finally, Sec. VII is 
devoted to our concluding remarks. 



II. THE TIME-DEPENDENT TMH 



Under the two-mode approximation, where the quantum field operators \l/ a 
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r, t) a 



36|, the 



and = <Pb{*,t)b are restricted to the fundamental states ipe (r, t) (£ = a,b) 
coupled Bose-Einstein condensates are described by the TD Hamiltonian (h = 1) 

H(t)=J2 ["t (t) + It (*) + lab if) atatfb 

£=a,b 

-g(t)(e~ iS ^a^b + e im atf), (1) 

where a and b are standard bosonic annihilatio n op erators, associated with condensation in 
hyperfine levels |2, 1) and |1, — 1), respectively [la . 
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37]] . The phase 5 (t) is associated to 



the detuning A(t) from the Raman resonance between the atomic transition |2, 1} <-> |1, —1), 
which may be a TD function (by varying the laser frequency), through the expression 5 (t) = 
A(r) dr+ 5q. The TD trap frequencies ue, the interspecies and intraspecies collision 



parameters jab and je, and the Rabi frequency g, follow from 
u t (t) = / d 3 rtp* e (r, t) 



^V 2 + V4(r,t) 



li(t) 



2m 



d 3 r\w (r,t)\ l 



(2a) 
(2b) 
(2c) 
(2d) 

where m is the atomic mass. We assume that the time dependence of the trap potential 
Vi (r, t) is generated by adiabatically varying the trapping magnetic field. Such adiabatic 



lab{t) = ^^ /d 3 r|^(r,t)^(r,t)| 2 



m 



9(t) = ^P- I dW a (r,t)^ b (v,t) 



variation of the trapping field has been assumed to ensure the validity of the two-mode 
approximation. The time- varying scattering lengths Anh(t) and Ae(t), are accomplished via 
Feshbach resonances, by tuning a bias magnetic field 38j. Finally, as mentioned above, in 
real experiments with atomic BECs the Rabi frequency may be a time- varying function since 
the amplitude and phase of the pumping fields may vary in time jl2j ]. 

Except for the Josephson-like coupling, the Fock states are eigenstates of all the terms 
in Hamiltonian ([I]). Thus, in order to get rid of this TD coupling in ([1]), we consider a 
transformation with the unitary operator 



V(t) = exp 



r(t) 



-i<t>{t) 



<)b) 



(3) 



analogous to that defined in Ref. 22]]) to obtain the transformed Hamiltonian 



H(t) = V^HV - iV^d t V = u&ht + K el (t) + H mel (t), 



(4) 



=a,b 



where ng = v-l is the number operator associated to each condensate having effective fre- 
quency 

u t (t) = uj e (t) + {25 lh -l)g (t) cos [0 (t) - 5 (t)} tan [r (t) /2] . (5) 

In the framework associated to the transformation ([3]), the system ends up with an inelastic 
collision term apart from the elastic one already present in (JTJ) . The Hamiltonians accounting 
for such interactions, also weighted by the TD function A(t) = -f a (t) + 7&(i) — J a b(t)i are 
given by 



n d {t) = | 7 a(t) cos 2 [r{t)/2] + 7fc(i) sin 2 [r{t)/2] 
+ ( 7o (t) sin 2 [r{t)/2] + lb {t) cos 2 [r(t)/2] 



A(t) 



sm 



V(t)]J(at)V 
r(t)]W)V 



sm 



+ { lab (t) + A(t) sin 2 [r(t)]} a)ab\ 
H mel (t) = sin[r( t)] - M sin[2r(t)] U-^) (a^ab 



+ < l7fc(t) ~ 7a(t)] sin [r(t)} + ^ sin [2r(t)]} e"^ aW 
+ M sin 2 [r(t)] ( e -^(t) a t 6 ) 2 + h .c. 



(6a) 



(6b) 
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The form of Hamiltonian (j4j) is established provided that the TD parameters r(t) and 4>(t) 
satisfy the coupled differential equations 



r(t) = 2g(t) sin [4> (t) - 5 (t)} , 

cj)(t) = uj(t) + 2g{t) cot [r{t)\ cos [0 (t) - 5 (t)} 



(7a) 
(7b) 



where 



to (t) = co a (t) - to b (t) . (8) 

The expression ([8]) represents an effective frequency for the system composed by the two- 
mode condensate, which plays an important role in the solutions of the characteristic equa- 
tions (J7j). In the Appendix we present a comprehensive analyzes of the analytical solutions 
of the Eqs. ([7]) for the on- and off-resonant regimes which are defined by comparing the 
effective frequency oj(t) with the detuning between the laser field and Raman transition 
A(t). The on-resonant regime, where A(t) = u (t), implies that the detuning from Raman 
transition must equals the effective frequency of the two-mode condensate. Otherwise, we 
have the off-resonant regime, where A(t) = oj(t) — to, w being some constant. 

Similar coupled differential equations were obtained by Smerzi et al. (39] in a semi- 
classical treatment of the double-well problem, and by Chen et al. j3] in a full quantum 
approach of the BECs in the two-mode approximation. In both references all the parameters 
in their Hamiltonians are constants, except for 5 (t) which, in Ref . .22] , is an adiabatically 
time- varying parameter. 



After the experiments by Hall and co-workers with 87 Rb [15|], where the scattering lengths 
satisfy the relation A a : A ab : A b = 1.03 : 1 : 0.97, and consequently 7 a ~ 7 b ~ 7^/2 (assum- 
ing spatial Gaussian function), a number of papers have driven attention to this particular 



case 



22 



401 ] whose Schrodinger equation is exactly soluble. However, if the proportion 



1.03 : 1 : 0.97 is broken, the system is not exactly integrable, but admits approximated solu- 
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22] . In our paper we are concerned with specific solutions of the 



tions as shown in Refs 

characteristic equations (J7J), under the rotating- wave approximation, which turn negligible 
the contribution of the inelastic interactions compared to the elastic one. This is done by 
substituting the solutions for r(t) and </>(t), obtained in the Appendix, into Eqs. §5§ and 
(Q, and rewriting Hamiltonian (j3J) in the interaction picture. Thus, after a time average of 
the TD parameters appearing in this Hamiltonian, we analyze the conditions leading to the 
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effective interaction 

Mr)n e + n d (r). (9) 

< ' t=a,b 

By adopting this procedure, where the inelastic interactions becomes despicable, we get he 
evolution operator U (t,t ) = exp ^— i f* 'Heff( T )d T ^ and, consequently, a prepared state 
(to)) evolves according to Hamiltonian ([T]) as 

|V (t)> =V(t)U (t, t ) (t ) IV (to)) • (10) 



III. DYNAMICS OF BECS FOR INITIAL BLOCH STATES 



Following Arecchi et al. [41[ and Dowling et a/. [42], we recall that the Bloch states (BS), 
also called atomic coherent states, spanned in the Dicke basis \j, m), where j = N/2 and 
\ m \ < j, are obtained through a specific rotation on the reference state 

\a,(3) = e^-°- i ^)\j,j), (11) 

where N is the number of condensed particles, J + , J_ (together with J z ), are the generators 
of the su(2) algebra, and (a, (3) are the polar and the azimuthal angles, respectively, defined 
on the Bloch sphere. The Heisenberg angular-momentum uncertainty relation for the BS 
reduces to 

((AJ x ,) 2 )((AJ y ,) 2 ) = \\(J z ,)\ 2 , (12) 

with the mean values being calculated in a rotated coordinate system x', y', z', where z' is 
an axis in the (a,/3) direction through the center of the Bloch sphere. Therefore, the Bloch 
vector is defined as the unit vector 

n = (sin a cos /3, sin a sin (3, cos a) , (13) 

in z'-axis. Using the Schwinger relations 

J x = g (^b + a6 f ) , (14a) 
J y = - {<Jb - atf) , (14b) 
J z = ~ (a f a - 6 f 6) . (14c) 

with J± = J x ± iJ y and the basis states {\j, m) = \N/2, (N a — Nb) /2)} = 
{|JV ) \N b ) = \N a , N b }}, where N a and A?" 6 (N = N a + N b ) stand for the number of atoms in 



the condensates, such that \j,j) •<=>- |iV, 0), it is straightforward to check that the BS can 
be defined through bosonic operators as 



\a,0) 



a 



cos I —j a) + sin [^] e l/3 b ] 



a 



N 



|0,0). 



(15) 



This state has a well-defined relative phase (3 between the two bosonic modes. 

Now, it is evident from relations ffl4|) that the unitary transformation V (t) turns to be 



exactly the rotation operator e~> 



-i0 



^ + \ if one considers r (t) = a and <j)(t) = (3. 



Therefore, it is easy to demonstrate through Eq. (flQj) that an initial BS \ip (to)) = |«o, A)) 
= |r o ,0o) evolves to another BS 



cos 



r(t) 



-iNtp N (t) 



-MvnW \ r (t),<f)(t)) 



a) + sin 



r(t) 



A' 



|0,0), 



(16) 



apart from the global phase factor e lN ^ N ^\ where 



<p N (t) = / {u a (r) + (N - 1) [ 7a (r) cos 2 [r(r)/2] 



to 



+ T&Msin 2 [r(r)/2] 



A(r) 



sin 2 [r(r)] 



dr. 



(17) 



The relative phase between the condensates 0(t) is associated to the mean values (J x (t)) = 
N sin [r(i)] cos [</>(t)] /2 and (J y (t)) = N sin [r(i)] sin [0(t)] /2, whereas r(t) is related to the 
population imbalance, AN(t) = (N a — N b ) = 2 (J z (t)), through the relation 

AN(t) = Ncos[r(t)\ . (18) 



We stress that according to the evolution operator, Eq. (flOl) . the evolved state (1161) 
remains a BS (apart from a global phase factor), since ao = r and (3q = 4>o- Note that 
the collision parameters are restricted to the global phase e~ tN '' PN ( t \ being irrelevant to the 
analyzes developed below of the nonadiabatic geometric phases acquired by the state vector. 
On the other hand, collisions become relevant when considering other initial states as the 
product of Glauber's coherent states \ao) \Po) instead of the BS 



IV. GEOMETRIC PHASES OF THE BS 



To study the evolution of the geometric phase in the two-mode BECs we use the kinematic 
approach developed by Mukunda and Simon |43| , where the geometric phase (pc is obtained 
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as the difference between the total phase 0r(£) = arg((0(to)| VK*))) an d the dynamical phase 
Mt) = ~i Jl (ip(r)\ £ |V»(r)> dr, resulting in 

9 



<V(r)| — |V(r))dr. 



(19) 



(j) G {t) = arg(^(t )| 

The expressions for an d <Ad(£) as function of the system parameters are, respectively, 

G (t) = jV a rg{cos (r /2) cos [r(t)/2] + e^ - *" 1 sin ( rQ / 2 ) s in [r(t)/2]} 

-- / 0(r) {1- cos [r(r)]}dr, (20) 



and 



<t* D {t) = -N [<p N (t) - <p N (t )] + 



N 



(t) {1 — cos [r(r)]} dr, 



(21) 



to 



where (pN^t) was given in Eq. (JT7I) . In Ref.|44j, the authors also obtain Eq. ([20l) . through 
the Gross-Pitaevskii equation, for the geometric phase acquired by the wave function of a 
BEC in the double-well problem, under the two-mode approximation. 

The time evolution of the BS can be followed on the Bloch sphere through the vector 
n (t) = (sin[r(t)] cos[0(£)], sin[r(t)] sin[0(i)], cos[r(t)]), for the different solutions r(t) and 
4>(t) presented in the Appendix. To follow such evolution and, consequently, to analyze its 
geometric phase 4>c(t), we must estimate the integrals in Eq.(j2]), i.e., the typical values for 
trap frequencies, Josephson-like coupling, intraspecies, and interspecies collision rates. For 
the sake of simplicity we model the effective frequency for both atomic species as harmonic 
trap potentials where the TD distribution for each condensate pi (r, t) can be approximated 
by a stationary Gaussian function such that 



1 



2irxl 



3/4 



-r 2 /4x^ 



(22) 



where Xt — ^/h/2muj£ stands for the position uncertainty in each harmonic oscillator ground 
231 ] . With this assumption the integrals in Eq.(j2J) are immedia tely estimated using 



state 



15 
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46|: m = 1.4 x 1(T 25 



typical physical parameters of the experiments with 87 Rb atoms 
Kg, ujt ~ 10 1-2 Hz, Ae ~ 5 nm, and fl ~ 10 3 Hz. To obtain some insights of the pseudo-spin 
dynamics under the TD Hamiltonian parameters, we consider the trap and Rabi frequencies 
as harmonic functions, oscillating around the typical constant values, as follow 



u e (t) = u eo + tot sin (xtt + &) 
git) = go + gsbn(fit), 



(23a) 
(23b) 
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with the parameters ueo, ue, xt-, 9o, 9, and /i being constant. Since the elastic collisions 
contribute only to a global phase factor, they will be assumed as the standard constant 
parameters in the literature 

7, = (24a) 

7- = — ■ (24b) 
m 

V. NONADIABATIC GEOMETRIC PHASES AND PSEUDO-SPIN DYNAMICS 
IN BECS 

In this section we present a detailed study of the geometric phase acquired by the whole 
wave function of the two-mode condensates. To this end, we plot the time evolution of the 
geometric phase f l20l) . and analyze its behavior through the evolution of Bloch vector (a 
map of the wave function of the BECs on Bloch sphere). This procedure allows for a better 
understanding and visualization of the concept of geometric phase for open trajectories in- 
troduced by Samuel and Bhandari [30| . In particular, we are interested in the dependence of 
the geometric phase on the constant of motions coming from the solutions of the character- 
istic equations (171) and also on the time dependence of the Hamiltonian parameters. In spite 
of the general solutions presented in the Appendix for these equations, we have assumed, 
for the analysis developed below, N = 1, 8$ = 0, and go = 6257T Hz. 



A. Solutions for r constant 



Before analyzing the geometric phases for the on- and off-resonant solutions, it is instruc- 
tive to present their evolutions for the simple case where the parameter r is kept constant 
while cf) (t) obeys Eq. fl37j) (since g — 0). In this case, the expression for the geometric phase 
coming from Eq. (I20l) . becomes 



Mt) = N jarg {cos 2 (r/2) + e -M'>-*»l sin 2 (r/2)} - iL^^l _ ^ j . ( 25 ) 

In Fig.l the absolute value for is plotted against the dimensionless r = u) a ot, 

assuming typical values u> a o = 2u>bo — Q2.5tc Hz. For r = tt/2, with the Bloch vector standing 
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on the equatorial plane, and u a = ujb = 0, the geometric phase evolves by jumps, as indicated 
by the thick solid line. These jumps occur every time the relative phase <ft connecting the final 
to the initial Bloch vectors equals (2n + 1) tt, n being an integer. The jump discontinuities 
occur because there are an infinite number of small geodesic-lengths connecting the vectors 
extremities, rendering the geometric phase undefined 48(. On the other hand, we observe 
that before jumping to <Pg{ t ) — i-e., for r < 2n, the geometric phase remain null since the 
small geodesic-length connecting the extremities equals the Bloch-vector trajectory itself. As 
soon as the Bloch vector acquires a relative phase larger than n, the small geodesic-length 
connecting the extremities completes a loop over the equator, making the acquired geometric 
phase proportional to nNn, where n, as defined above, turns out to be the winding number, 
i.e., the number of loops around the z-axis of the sphere. 

The same interpretation given above for the geometric phase holds for the two other 
curves obtained for r = vr/2.1, except that the jump discontinuities are substituted by high- 
slope curves around the points where = (2n + 1)tt. Moreover, the net effect coming from 
the TD parameters of the Hamiltonian is to delay or advance the sequential increments 
of the relative phase <fi and, consequently, of the geometric phase, as observed from Fig.l. 
The solid line, obtained for r = tt/2.1 and u a = uj^ = 0, shows that the increments of the 
geometric phase, besides being smaller, present the same rate of variation when compared to 
the case r = tt/2. When the trap frequencies are oscillating functions, with u a = u h = u a0 /4, 
Xa = Xb = £ a = 0, and = 7r/2, the time- dependence shown by the dashed line, 

results in the oscillations of the time intervals between the increments of the geometric phase. 

To better visualize the above discussion about geometric phases acquired in open tra- 
jectories, in Fig. 2 we plot the evolution of the Bloch vectors for the cases r = tt/2 and 
r = 7r/2.1, with u a = cj& = 0, considering the same time interval r = Air. The black and 
grey vectors indicate the coincident positions of the initial and final Bloch vectors, for the 
cases r = 7r/2.1 and r = tt/2, respectively, after a complete rotation around the sphere 
whose directions are indicated by the arrows. The trajectories described by the black and 
grey vectors are indicated by solid and dashed curves, respectively. Evidently, the geometric 
phase acquired during the evolution of the Bloch vector in the case r = tt/2.1 (the solid 
angle comprehended by the semi-hemisphere above the solid circumference), is smaller than 
that for the case r = tt/2 (the solid angle corresponding to the north hemisphere, equal to 
2tt). 



11 



The solution with r constant means steady population imbalance AN, whereas the rel- 
ative phase (pit), another parameter examined by experimentalists and necessary to com- 
pletely define the BS, is a linear function of time when u a = ujf, = or an oscillating function 
when u a = ujfy = u> a o/4. Note that the dynamics of the population imbalance and the relative 
phase may be followed through the projection of the Bloch vector trajectory on z-axis and 
x-y plane, respectively. 



B. On-resonant solution 



As indicated in the Appendix, through the constant of motion C = 
sin [r(t)} cos [4>(t) — 5(t)} we obtain the solution of the characteristic equations ([7]) in 
the on-resonant regime where A(t) = uj{t). All the possible trajectories in the portrait 
space r(t) x (4>(t) — S(t)), are restrained to the level curves obtained as projection of the 
surface plotted in Fig. 3, which follows from C. 

To better understand the geometric phase acquired by the state vector |\l/(t)) in the 
on-resonant solutions, we consider two different cases, A = and A ^ 0, and analyze its 
dependence on the constant C, through the relation 



4> G (t) = iVarg {cos (r /2) cos [r (t) /2] + e^^ 01 sin (r /2) sin [r (t) /2]} 

2Cg(t') cos [r(t')\ 



to 



1. The case A = 



^Jdt'S [ u;(t'){l-cos[r(t')}} + ^ 



+ cos [r (£')] 



(26) 



In Fig. 4 we plot the evolution of the geometric phase against r = g t, considering uj a = 
ujb = g = 0. The thick solid line on the abscissa axis corresponds to the choice r = vr 
and 0o = 77 /2, leading to C = 0, under which the geometric phase is null or undefined 
as indicated by the open dots over the abscissa axis. Note that for r = vr we get, at 
t = 0, an undetermined equation ( I31b[) for <p{t). To circumvent such indetermination we 
impose on Eq. (1301) the constraint <ft(t) — 5(t) = (2n + l)7r/2 over any time interval, to 
determine 4>(t) independently of Eq. (131bl) . Since for A = it follows that 5(t) = 5 , 
implying that <p{t) = 5 + (2n+ 1)tt/2, the geometric phase for the case C = A = simplifies 
to (pcif) — N aig {cos [(r it) — r ) /2]} and, consequently, (pcif) is nuu f° r l r if) — r o| < n 
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and undefined for \r(t) — r | = n. Still in Fig. 4, the solid and dashed lines, associated 
to the pairs (r , 4>o) = (tt/5, 7r/4) and (ir/A, 37r/10), respectively, correspond to the same 
constant C ~ 0.41. These curves exhibits similar behaviors due to the fact that, with the 
same constant C, they present the same trajectory on the portrait space of Fig. 3, despite 
starting from different initial conditions. The dotted and dashed-dotted lines, associated to 
the pairs (vr/3, 0) and (vr/3, n), and corresponding to the constants C ~ 0.87 and C ~ —0.87, 
respectively, are symmetric around the abscissa axis r. Such property of symmetry reflection 
of the geometric phase (fie — (fie, is consequence of the change O ~ ¥ 0o ± tt, implying that 
C — > — C. It is worth noting that the larger the absolute value of C the smaller the acquired 
geometric phase and vice-versa. 

In Fig. 5 we plot the evolution of the Bloch vectors coming from the on-resonant solution 
with the initial conditions (ir, 7r/2) and (n/3, 0) corresponding to C = and C ~ 0.87, 
whose initial and final positions are represented by the black and grey vectors, respectively. 
As in Fig. 2, we consider the evolution of both vectors during the same time interval r = n. 
Through the solid line trajectory described by the black vector, which oscillates between 
the north and south poles, it is straightforward to conclude that the geometric phase is null 
during the whole time evolution, except when the vector reaches the north pole, where the 
geometric phase becomes undetermined. As the dashed trajectory of the grey vector is not 
restricted to a meridian, as in the case C = 0, the geometric phase acquired is evidently 
non-null. 

2. The case A / 

As we are interested in the dependence of the geometric phase on constant C and, now, 
on the effective frequency of the two- mode condensate A(t) = cu(t), we consider all the 
parameters of the Hamiltonian being time- independent, u a = ^fe = g = 0, except for 
S(t) = S + uj{t')dt' . In Fig. 6 we plot the geometric phase against r = g t for different 

J t 

initial conditions (r , 0o) an d effective frequencies A. As indicated by the thick solid line 
associated to the initial conditions (ir, vr/2), corresponding to C — 0, with u a = 2u b = go/10, 
the geometric phase is not null, differently from the case A = 0. The discontinuity exhibited 
by this curve follows from the arg function, whose characteristic jumps occurs whenever 
r n = [{2n — 1) 7r<7o ] /2cu, n being a positive integer. As indicated by the solid and dashed 
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lines, the property of symmetry reflection of the geometric phase in the abscissa axis still 
follows when changing, simultaneously, (po — > 4>o ± 7T, and to — cj. In fact, the solid 
line corresponds to the initial condition (vr/3, 0) with C ~ 0.87 and u a = 2u>b = go/10, 
while the dashed line corresponds to (tt/3, it), with C ~ —0.87 and uj, = 2u a = go/10. 
When the change 0o — > 0o =t ^ are n °t followed by that a; — > — a;, such a symmetry is not 
accomplished as indicated by the dotted line corresponding to the initial conditions (ir/3, n) 
with C ~ —0.87 and, now, u a = 2u>b = go/10. 

To visualize the acquisition of the geometric phase we again return to the evolution of 
the Bloch vector. As observed from Fig. 7, the solid trajectory described by the vector 
associated to C = 0, whose coincident initial and final positions are indicated by the black 
vector, leads to a finite solid angle and, consequently, a finite geometric phase (during the 
time evolution r = n considered for both cases presented). The control of this solid angle 
may be accomplished through the parameter A — the larger A the larger the solid angle and 
vice-versa — as demonstrated experimentally through the polarization vector of a photon 
undergoing a Mach-Zehder interferometer [47J . The evolution of the grey vector, associated 
to C ~ 0.87, exhibits a dashed trajectory where the initial and final positions are slightly 
different. As the time evolution proceeds, such trajectory leads to a geometric phase which 
easily exceed that of the case C = 0, as indicated in Fig. 6. 

The population imbalance AN for the cases A = and A / is a time-oscillating 
function strictly dependent on the shape of Rabi frequency g(t), thus exhibiting a strong 
connection between the dynamics of population inversion in two-level systems and popula- 
tion imbalance of BECs. In fact, when written the Hamiltonian (TjQ) through the quase-spin 
operators in Eqs. fll3l) . we obtain, apart from the collision terms, a driven interaction. On 
the other hand, the relative phase <p(t) depends also on the detuning S(t) besides g(t), as 
shown by Eqs. lpTT) . 

C. Off-resonant solution 

To analyze the off-resonant solution, where the detuning A(t) = oj(t) — to is controlled by 
adjusting the parameter w, we impose a constant Rabi frequency g(t) = go which implies a 
constant of motion C = rj sin [r(t)] cos [<p(t) — 5(t)] — cos [r(t)}. Similarly to the on-resonant 
case, all the possible trajectories for the off-resonant r(t) and <p(t) are restrained to the 
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level curves of the surface following from C, presented in Fig. 8. When 77 = 2g /zu 3> 1 it 
is verified that the constant C reduces to that of the on-resonant solutions, unless for the 
multiplicative factor 77, and the surface presented in Fig. 8 also reduces to that of Fig. 3. 
However, for rj <C 1 we obtain an approximated constant value of r. Finally, when rj ~ 1, we 
obtain the surface whose level curves encapsulate all the possible trajectories in the portrait 
space r(t) x <f>(t) — 5(t), as shown in Fig.8. 

To study the effect of rj on the off-resonant geometric phase, given by 

<f> G (t) = iVarg {cos (r /2) cos [r (t) /2] + e ~ i[<Kt) - M sin (r /2) sin [r (t) /2]} 
t 

JV / , , f . r , ,,,, r , ... C + cos \r (t')} 

- -J it [.(f) {1 - cos [r (*<)]} + ro eos [r (0] 1+cos| ;,^ 

to 

we plot in Fig.9 the evolution of for different values of 77, all starting from the point 

(7r/3, 0), with (j a o = 2ubo = <?o/10 and cD a = ujf, = 0. The dotted line, following from the 
case where 77 = 40, indicates a similar behavior to the corresponding case of the on-resonant 
solution (represented by the solid line in Fig. 6). For rj = 0.1, the solid line shows that 
the geometric phase is a strongly oscillating function, a behavior that is better visualized 
through the evolution of the corresponding vector in Bloch sphere, Fig. 10. Finally, when 
zu ~ g , as for 77 = 2, the geometric phase shows discontinuities, as indicated by the thick 
solid line, which turns out to be a signature of the off-resonant solution. We note that the 
property of reflection exhibited by the geometric phase coming from the on-resonant solution 
is also present here when substituting, simultaneously, 0o ~~ ^ <Po ± tt, uo — > — u>, and 77 — > —77. 

In Fig. 10 we present the evolution of the Bloch vectors in the time interval r ~ 97r/10, 
for the cases 77 = 0.1 and 77 = 2, both starting from the common point (tt/3, 0). The black 
vector for 77 = 0.1, presents a behavior limited to the north hemisphere, described by the 
solid trajectory, which exhibits periodic up and down motions on the parallels, responsible 
for the oscillations of the geometric phase showed in Fig.9. The grey vector for 77 = 2, 
by its turn, indicates a rather complicated dashed trajectory which descends to the south 
hemisphere and goes back to the north. 

Besides depending on the Rabi frequency, as in the on-resonant case, the population 
imbalance for the off-resonant solution also depends on the detuning zu which makes its 
mean value not null. The relative phase depends on S(t), g(t), and the detuning zu. 
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D. Time-Dependent Effects on the Geometric Phases 

The behavior of the geometric phase when both, the trap and Rabi frequencies are TD 
harmonic functions, as described by Eq. (1231) . is analyzed in Fig. 11, where we plot <Pg{ t ) 
against r = g$t, considering the same initial conditions (vr/3, 0) for all the curves. Starting 
with the resonant solution with A = 0, we obtain the dotted curve for the parameters 
u a = u b = and g = fi = g , to be compared with the dotted curve of Fig.4. We observe that 
both dotted curves are very close to each other, with the the increasing rate of the geometric 
phase being modulated by the TD Rabi frequency in Fig. 11. The solid line corresponds to 
the on-resonant solution with A ^ 0, for the parameters u a0 = 2u b0 = go/10, co a = u b = 0, 
and g = fi = go. This curve is to be compared with the solid line in Fig. 6, showing again 
that the increasing rate of ipc can be controlled through the TD Rabi frequency. The dashed 
line, also corresponding to A ^ 0, with u a0 = 2u b0 = g , u a = u b = 0, and g = /i = g , 
shows that the increasing rate of the geometric phase may also be controlled through the 
trap frequencies. Finally, the thick solid line, corresponding to the off-resonant solution, 
with uj a0 = 2u m = g , u> a = u b = \a = Xb = go/2, £a = 0, £ b = -ir/2, and g = 0, to 
be compared with the thick solid line of Fig.9, also indicates the important role played by 
the time dependence of the trap frequency in the geometric phase. We finally observe that, 
evidently, these TD effects have direct implications on the population imbalance and relative 
phase between the condensates. 



VI. GEOMETRIC PHASES AND THE EXTERNAL JOSEPHSON EFFECT 



In this section we analyze the two-mode condensates from a different perspective: as 
a single atomic specie trapped in a symmetric or asymmetric double-well potential. As 
the internal Josephson effect is here substituted by the tunneling interaction, the laser 
pumping becomes unnecessary and we impose that A = 0, such that 5 = 0. Moreover, 
in the external Josephson effect the interspecies collision rate correspond to a second order 
correction compared to the intraspecies collision rates, justifying the assumption that ^ab — 



36|,|49|. 



The Hamiltonian ([T]) applied to this different physical situation leads, under the above 
restrictions, to a similar transformed interaction (j4j) and characteristic equations (JTj). There- 
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fore, the different solutions of the coupled differential equations apply directly to the external 
Josephson effect, with the on- and off-resonant processes describing the symmetric (u = 0) 
and asymmetric (a? 7^ 0) wells solutions, respectively. 



A. Symmetric wells 

^From the above discussion we readily verify that the solutions for r(t) and </>(£), com- 
ing from the symmetric wells, are give by Eqs. (1311) and ( 1321) with the constant of motion 
C = sin [r(t)] cos [</>(£)]. We observe that the TD tunneling rate g(t) is accomplished by 
modulating the amplitude of the counter propagating classical fields that generate the bar- 
rier. Similarly to the internal Josephson effect, all the possible trajectories for r(t) and 4>(t) 
follow from the level curves of the surface plotted in Fig. 8, assuming 5 = 0. The same level 
curves were obtained, by numerical methods, in Ref . [44j . The geometric phases acquired 
by the evolution of the state vector of the BECs are given by Fig. 4 and the Bloch-vector 
trajectories by Fig. 5, both obtained from the on-resonant solution of the internal Josephson 
effect with 5 = 0. 



B. Asymmetric wells 

The solutions of the characteristic equations for the asymmetric wells follow by imposing 
constant values for the effective frequency u = w (since A = 0) and the tunneling rate 
g = go. We thus obtain the solutions (1341) with rj replaced by rj = 2g /u and 5 = 0. The 
constant of motion becomes C = rj cos [4>(t)] sin [r(t)\ — cos [r(i)]. The phase-space portrait 
r(t) x cj)(t) is given by Fig.8 (with 5 = 0), whose level curves indicate all the possible 
trajectories for r(t) and 4>(t). The same level curves were obtained by numerical methods in 



Ref. [44j . As an example of the geometric phase acquired by the evolution of the state vector 
in the asymmetric wells, we take the dotted line curve of Fig. 9, corresponding to the case 
uj = w = go/20. (The other two curves in Fig.9 do not satisfy the condition u = w.) The 
trajectory of the Bloch vector associated to uj = go/ 20 is approximately give by the dashed 
curve in Fig. 7. 
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VII. CONCLUDING REMARKS 



In the present work we analyze the dynamics of two interacting condensates, with a full 
TD Hamiltonian. Starting from the Hamiltonian ([T]) under the two-mode approximation, 
an effective interaction Q is established under the RWA, provided that the polar r(t) and 
azimuthal <p(t) angles, which define a Bloch state, satisfy coupled differential equations. 
This procedure enable us to define a detuning A(t) from the Raman resonance between the 
atomic transition, together with an effective frequency for the condensates u(t), as in Eq. 
(jSJ). Thus, two different solutions arise for the differential equations coupling the parameters 
r(t) and <f>(t), the on-resonant solution, where A(£) = u)(t), and the off-resonant solution 
where A(i) = to (t) + m, w being a constant. After solving analytically the coupled equations 
for both regimes, we present a detailed analyzes of the geometric phases acquired by the 
Bloch state of the system, also discussing the relative phase and population imbalance. 

A main result of our work is the connection between geometric phases and constant of 
motions of the interacting condensates which are identified through the analytical solutions 
of the coupled differential equations for r(t) and <f>(t). For each on- or off- resonant solution 
we assign a constant of motion which determine the dynamical behavior of the state of the 
system and, consequently, its geometric phase. We also note that these constants of motions 
follow from level curves in the portrait space which are also obtained analytically. 

To better visualize the time-evolution of the geometric and relative phases acquired by the 
state vector, together with population imbalance, we also analyze the trajectory of the state 
vector mapped on Bloch sphere. Finally, we present a brief discussion of the TD effects of 
the trap and Rabi frequencies in the geometric phases together with the connection between 
its evolution in both cases of internal and external Josephson coupling. 

As in this work we studied only the evolution of an initial Bloch state, the collision pa- 
rameters were restricted to the global phase of the evolved state, which remains a Bloch state 
under the two-mode approximation. Therefore, collisions terms, which are also assumed as 
TD parameters, do not play a decisive role in our present analyzes. It is worth to consider 
distinct initial state to analyze the effects of the collisions parameters in the dynamics of 
the geometric phase and population imbalance. 
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VIII. APPENDIX 



A. Analytical solutions of the characteristic equations ([7]) 

In this Appendix we present some specific solutions of the characteristic equations (j7j), 
following a more detailed treatment in 50j. We investigate two different regimes of the laser 
field amplification, the on-resonant and off-resonant regimes, which are defined comparing 
the effective frequency of the two-mode condensate, u(t), with the detuning between the 
laser field and the Raman transition A(£). As mentioned above, in the on-resonant regime, 
where A(£) = u (t), the rate of time variation of the laser field equals the effective frequency 
of the two-mode condensate. Otherwise, we have the off-resonant regime. 

B. On-resonant process 

Defining x(t) = </>{t) — S(t), the characteristic equations (JTj) becomes 



r(t) = 2g(t) sin [ X {t)\ , 

X (t) = u(t) - A(t) + 2g(t) cos [ X {t)] cot [r(t)] 



(28a) 
(28b) 



such that, in the on- resonant regime we are left with the first-order differential equation 

(29) 

After integrating Eq. (129"]) we obtain the constant of motion 



dr 

— = tan x tan r. 

d X 



sin [r(t)} cos [<j>(t) - 5(t)] = C, 



(30) 



with C depending on the initial values r , O , and 5 . Thus, the resonant solutions of 
Eqs.f T2"8"j) . are given by 

cos r 



cos [r(t)\ = Vl — C 2 si 



sin 



u(t,t Q ) + arcsin 



(p(t) = 5(t) + arccos < — 



c 



where 



sin [r(t)\ J ' 



u(t,t ) = -2 / g(r)dr. 

J to 



(31a) 
(31b) 

(32) 
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C. Off-resonant process 



Considering the off-resonant regime, where A(£) = u(t) — zu, w being a constant, Eqs. 
can again be solved by quadrature as far as we assume the Rabi frequency g to be also 
a constant g . In this regime, defining 77 = 2g /zu, Eq. (l3"0l is replaced by 

77 cos {(j){t) - 5(t)] sin [r(t)] - cos [r(t)] = C, (33) 

which again depends on the initial values r , O , S . The solutions for this case are given by 



cos [r(t)} 



sin 



+ arcsm 



1 +V 2 

1 + rf) cosro + C 



+ 7] 2 (t- t ) 

c 



(f)(t) = 5(t) + arccos 



Tjy^l + if — C 2 

C + cos [r(t)) 



I + 77 2 



?7sin [r(t)} 



(34a) 
(34b) 



D. A constant solution for r 

Another solution for r(t) and <p(t) arises when we impose that r remains constant in time. 
Through this solution, given by 

r(t) = r , (35) 

0(t) = S(t) + rrn = 0o + / Hr) + 2(-) n g (r) cot {r )]dr, (36) 

with r 7^ n7r and n being an integer, the state vector of the system acquire only relative 
phase 0(t). This formal solution for 0(t) can be even simplified noting that the physical 
implementation of this regime requires, necessarily, that g — 0. In fact, for the population 
imbalance to be null, due to the constant value for r, the Rabi frequency must also be null, 
therefore, Eq. fl36|) simplifies to 

0(i) = 0o + / co(r)dr. (37) 

J t 
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Figure captions 

Fig. 1 Absolute value of geometric phase (pcif) against the dimensionless r = uj a ot, for r 
constant, with u> a0 = 2u b0 = 62.57T Hz. 

Fig. 2 Evolution of the Bloch vectors coming from constant solutions of r in the time 
interval r = 4vr. The grey and black vectors correspond to the initial conditions (jr/2,0) 
and (7r/2.1, 0), respectively, with u a o = 2ubo = Q2.5n Hz and u a = u b = 0. 

Fig. 3 The portrait space r(t) x <p(t) — 5(t) obtained as projection of the surface which 
follows from the on-resonant constant of motion C = sin [r(t)] cos [4>(t) — S(t)]. 

Fig. 4 Evolution of the geometric phase (pcif) against r = g t, for on-resonant solutions 
of the characteristic equations (171) , with A = and u a = uj b = g = 0. 

Fig. 5 Evolution of the Bloch vectors coming from the on- resonant solutions of Eqs.([7|) 
in the time interval r = 71. The black and grey vectors correspond to the initial conditions 
(•7T, 7r/2) and (vr/3, 0), with A = and u a = u h = g = 0. 

Fig. 6 Evolution of the geometric phase (pcif) against r = got, for on-resonant solutions 
of the characteristic equations (GO), with A ^ and co a = ujb = g = 0. 

Fig. 7 Evolution of the Bloch vectors coming from the on- resonant solutions of Eqs.([7j) 
in the time interval r = 71. The black and grey vectors correspond to the initial conditions 
(n, 7r/2) and (vr/3, 0), with A ^ and u a = £>& = g = 0. 

Fig. 8 The portrait space r(t) x <p(t) — 5(t) obtained as projection of the surface which 
follows from the off-resonant constant of motion C = 77 sin [r(t)} cos [4>{t) — S(t)] — cos [r(t)}. 

Fig. 9 Evolution of geometric phase 0g(£) against r = g t, for off-resonant solutions of 
the characteristic equations (J7J) and different values of r], all starting from the point (tt/3, 
0), with u a0 = 2u b0 = g / 10 and u a = u b = 0. 

Fig. 10 Evolution of the Bloch vectors coming from the off-resonant solutions of Eqs.([7|) 
in the time interval r ~ 97r/10. The black and grey vectors, corresponding to rj = 0.1 and 
r] = 2, respectively, both start from the common point (vr/3, 0), with u a0 = 2uJbo = go/lO 
and u a = ujb = 0. 

Fig. 11 Evolution of geometric phase <Pg{ t ) against r = got, for on- and off-resonant 
solutions of the characteristic equations ([7]), considering the same initial conditions (vr/3, 0) 
for all the curves. 
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